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^ . Abstract 

We investigate M — A SYM coupled to fundamental flavours at nonzero imaginary 
►^ , quark chemical potential in the strong coupling and large iV limit, using gauge/gravity 

t^ . duality applied to the D3-D7 system, treating flavours in the probe approximation. The 

'NT ■ 

0^ ■ interplay between Zjv symmetry and the imaginary chemical potential yields a series 



of first-order Roberge- Weiss transitions. An additional thermal transition separates 

in ■ , . . . 

f^ ' phases where quarks are bound/unbound into mesons. This results in a set of Roberge- 

Weiss endpoints: we establish that these are triple points, determine the Roberge- Weiss 
temperature, give the curvature of the phase boundaries and confirm that the theory is 
analytic in /i^ when fi^ w 0. 



1 Introduction 

The duality between strongly interacting gauge theories and gravity, via the AdS/CFT cor- 
respondence [1], has fuelled the exploration of a large number of strongly coupled systems 
and their associated thermodynamic phase structures. Essential features of the thermody- 
namics of pure gauge theories such as the deconfinement transition at strong coupling are 
captured by semiclassical gravitational transitions [2,3]. When applied to AA = 4 supersym- 
metric (SUSY) gauge theory, the gauge/gravity correspondence yields the full phase diagram 
at strong coupling, as a function of both temperature and chemical potentials for global 
R-symmetries [4-7]. Various qualitative aspects of this phase diagram are also mirrored by 
the weakly coupled gauge theory [8-12]. Extensions of such holographic computations to 
include gauge theories with matter fields transforming in the fundamental representation of 
the gauge group have been considered [13, 14] and the resulting phase diagrams have been 
extensively studied by various groups [15-28]. 

In this paper we study the thermodynamics of the SU{N) J\f = 4 theory coupled to 
M = 2 SUSY matter with an imaginary chemical potential for quark (flavour) number, 
using the D3-D7 setup at strong coupling, introduced in [13]. The theory has Nf flavour 
multiplets preserving M = 2 SUSY, with Nj <^ N, so that flavour loops are suppressed 
in the large- A^ theory. The physical motivation and rationale behind the investigation of 
thermodynamics with an imaginary chemical potential fij is discussed below; in short, it 
provides an additional window into the dynamics and phase structure of "QCD-like" gauge 
theories with fundamental matter fields. 

Our study of the gravity dual description of the strongly coupled gauge theory reveals the 
phase diagram depicted in Figure 1. At high temperatures T, there is an infinite sequence 
of first order transitions at m = /Urw > where 

ttT 
/XRW = (2r - 1)— , rGZ, (1) 

characterized by jumps in the phase of the Polyakov loop order parameter. The quarks, 
i.e. the fundamental matter degrees of freedom, remain unbound and are represented by 
the so-called black-hole embeddings of probe D7-branes in the gravity dual setup. The 
sequence of first order transition lines appears due to the interplay between Z^v-invariance 
of the adjoint fields and the imaginary quark chemical potential m, as was argued on 
general grounds by Roberge and Weiss in the context of QCD [29]. It reflects the presence 
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Figure 1: Phase diagram in the plane of temperature T and imaginary quark chemical potential 
/i/. The blue lines represent the first order Roberge- Weiss transitions in the high-temperature 
phase. The red line is the line of first-order phase transitions separating the high- and the low- 
temperature phases. The Roberge- Weiss endpoints, indicated by the black squares, are triple points 
at ifJ^mv , Tyiw) ■ The phase of the Polyakov loop is denoted with a. 



of unbound or deconfined fundamental matter degrees of freedom (carrying non-zero A^- 
ality) . A significant feature of the Roberge- Weiss (RW) lines is that they must end at some 
low temperature, below which quarks are confined or bound, thus acting as a diagnostic 
for confinement/deconfinement of quarks. The order of the transition at the Roberge- 
Weiss endpoints depends on the theory in question. In the holographic example we study 
here, the endpoints are triple points since the RW lines meet another line of first order 
transitions (indicated in red in Figure 1) at T = Trw. Below this line the quarks are 
bound into mesons corresponding to fluctuations of the "Minkowski-embeddings" of probe 
D7-branes in the gravity dual [30]. The line of meson-melting transitions is an extension 
to imaginary values of the same phenomenon found in this model for real values of the 
chemical potential [15, 18,21,26]. We find that this phase boundary is given by a set of 
parabola, centered at m/T = Inr/N (r G Z). Explicitly, we find 

—J , ^«43.5, (2) 

when 

(2r- - l)n LLi (2r + l)7r 

AT < T < AT • ^^^ 

Here Tq = Tc{iii = 0) is the critical temperature for vanishing /i/ and A = gyj^^fN is the 

't Hooft coupling, which is taken to be large. The Roberge- Weiss endpoints are therefore 

located at (/iRWi^Rw)) where /^rw is given by Eq.(l) and 
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The meson melting temperature at zero chemical potential is determined by the typical 
meson mass scale M in this model and Tq ~ 0.77M, where M = 2mq/v\ with vriq the 
common mass of the N = 2 flavour multiplets [26,30]. 

Since the flavour degrees of freedom are probes in the large- iV theory, they do not affect 
the dynamics of the adjoint degrees of freedom of AA = 4 SYM, which remain deconfined 
at all temperatures. Hence, the low-temperature phase does not resemble the gluon sector 
of QCD. However, fundamental matter in this system is bound into mesonic states at low 
temperature and unbound at high temperature (with the scale set by the meson mass). A 
notion of a low- and a high-temperature phase therefore still exists. 

The primary motivation for the exploration of thermodynamics in the presence of imag- 
inary chemical potentials comes from QCD. As is well known, the phase structure of QCD 
at nonzero temperature and quark chemical potential is difficult to determine. At large 
temperature and/or chemical potential predictions can be made with perturbation theory, 
due to asymptotic freedom, but most of the phase diagram requires a nonperturbative eval- 
uation. Thermal properties (T > 0) at vanishing chemical potential (/U = 0) are being 
investigated vigorously using lattice QCD; for a recent review, see e.g. [32]. At non-zero 
chemical potential however, a straightforward application of lattice QCD runs into prob- 
lems. Since the fermion determinant is complex, algorithms based on importance sampling 
can no longer be applied, resulting in the infamous sign problem. Various approaches to 
circumvent the sign problem in QCD at finite density have been developed over the past 
decade; an up-to-date review can be found in [33]. 

The sign problem is due to the nonhermiticity of the Dirac operator. Let Ip denote the 
massless Dirac operator, satisfying ^^Ip^^ = Ip- Integrating out the fermion fields then 
yields the determinant det M{p) = det(^ + m + /U70). Using 75-hermiticity it immediately 
follows that this determinant is complex, 

[det M(//)]* = det M{-fi*), (5) 

unless the chemical potential vanishes or is purely imaginary. This last observation suggests 
that lattice QCD simulations at imaginary chemical potential n = i^i can be performed 
using standard algorithms. Information about the phase structure at real chemical potential 
can then be obtained by analytic continuation: since the partition function is an even 
function of /U, analytical continuation from imaginary ^ [fi^ < 0) to real /U {fi"^ > 0) is in 
principle straightforward. This approach has indeed been used extensively [34-43]. 



As mentioned above, the phase structure at nonzero imaginary chemical potential is 
nontrivial due to the interplay between the Zjv symmetry and fij [29]. As we will review 
below, the partition function is periodic in nj and satisfies 

27rrT" 



Z[fii] = Z 



/U/ + 



N 



r£Z. (6) 



Combined with reflection symmetry Z[fii] = Z[—fii], this leads to Roberge- Weiss transitions 
at /[// = /Urw (Eq.(l)). At these values of fij, the theory undergoes a transition between 
adjacent Zn sectors, which can be distinguished by the phase of the Polyakov loop. 

At high temperature the Roberge- Weiss transitions are first order, resulting in non- 
analytic behaviour of observables at fij = /Urw- At low temperature, the Roberge- Weiss 
periodicity is smoothly realized and no nonanalyticity is present. This was first suggested 
in [29] and has been confirmed with lattice QCD simulations [34-43]. The Roberge- Weiss 
lines therefore end. In QCD the high- and low-temperature phase are separated by a thermal 
confinement/deconfinement transition,^ which extends into the fij direction. It is widely 
assumed that this transition line will join the Roberge- Weiss lines at m = //rw- Generically 
the phase diagram will therefore be as depicted in Fig. 1. 

In QCD properties of the confinement/deconfinement transition depend crucially on the 
number of fiavours and the quark masses. For instance, for three degenerate flavours, the 
transition is first order for light and heavy quarks (which is tied to the chiral and centre 
symmetries respectively), while for intermediate quark masses, the transition is a crossover. 
It has only recently been demonstrated that the Roberge- Weiss endpoint is a triple point 
for two [42] and three [37] degenerate light or heavy fiavours, while it is a second order 
endpoint in the 3d Ising universality class for intermediate quark masses. Moreover, in a 
beautiful paper [37], de Forcrand and Philipsen argued in the case of heavy quarks that 
the tricritical point found in the quark mass-temperature plane at fij = ^rw dictates the 
deconfinement critical line for real fi by tricritical scaling. 

It is therefore clear that an understanding of the phase diagram at imaginary chemical 
potential can provide important insight into the phase structure at real chemical potential. 
For this reason this topic has been under intense investigation, using both lattice QCD as 
mentioned above, renormalization-group studies [44], as well as effective models [45-48]. In 
this paper we take a third approach and apply gauge-gravity duality to study the phase 
diagram of a strongly coupled gauge theory at imaginary chemical potential. 



'^ Although strictly speaking it is only a confinement/deconfinenient transition when the quarks are in- 
finitely heavy, we use this terminology for fight quarks as weU. 



This paper is organized as follows. We begin by reviewing the arguments of [29] in 
Section 2 and summarize results that can be obtained at arbitrarily weak gauge coupling. 
In Section 3 we review the well known properties of the gauge theory model we consider 
(M = 4 SYM coupled to AA = 2 matter). Section 4 contains the holographic description of 
the model at imaginary chemical potential and a numerical study of the configurations that 
dominate the ensemble as a function of temperature and ^u/. The resulting phase diagram 
is constructed in Section 5. Section 6 contains a brief outlook. In the Appendix we outline, 
in some detail, the thermodynamics of the model at weak coupling. 

2 Roberge- Weiss phase transitions 

It is well known that the deconfinement phase transition in SU{N) gauge theories without 
matter in the fundamental representation can be precisely characterized, in the Euclidean 
formulation, by the expectation value of the thermal Wilson line or Polyakov loop P = 
jjTrexp{i § Ar). The Tjn symmetry of the pure gauge theory under which P — )■ p^'^'^^'^/" ^ 
{k = 0, 1, 2, ... A^ — 1) is spontaneously broken by an expectation value for the Polyakov 
loop in the high temperature deconfined phase. This distinction between confining and 
deconfined phases is lost upon the introduction of matter in the fundamental representation 
of the gauge group since there is no Zjv symmetry and the operator P can obtain an 
expectation value at all temperatures. However, as was argued in [29], there is a remnant 
of this symmetry which can still be used to distinguish between the phases in which quarks 
are confined or deconfined. 

The idea of [29] is that in theories with fundamental matter and an appropriately defined 
baryon number U{\)b symmetry, the difference between the confined and deconfined phases 
can manifest itself in the behaviour of the theory as a function of an imaginary chemical 
potential for baryon or quark number. With an imaginary chemical potential /x/ for quark 
number, the thermal partition function 

Z[^i] = Tr L-PH+ipt^iNA ^7^ 

is naturally periodic under m ^ m + 2ttT (/3 = 1/T), since Nq is quantized. Hence 
f G[-vr,7r]. 

There is a further periodicity of this partition function following from the Zjv symmetry 
of the adjoint degrees of freedom [29] and it can be understood as follows. We recall that the 
(imaginary) chemical potential for fermion number couples to fundamental matter fields in 
the same way as the time component Ar of the SU{N) gauge field. The imaginary chemical 



potential /i/ is introduced by the replacement dr + iAr — )■ 5,- + iA-r — i^j , acting on matter 
fields in the fundamental representation. 

This ;U/-dependence at the level of the action can be removed by a phase redefinition of 
fields (f) charged under U{1)b-, and transferred into boundary conditions around the thermal 
circle, 

<A(f,/?) = ±e^^^/^0(x,O), (8) 

where the + and — signs apply to bosonic and fermionic modes respectively. We can perform 
a second variable change, corresponding to a gauge transformation by an element U{x, r) 
of SU{N), with the property U{x,f3) = e^'^*^'^ [/(x, 0), r G Z. This leaves the action and 
the path integral measure invariant, whilst only altering the boundary conditions for the 
fields around the thermal circle, 

^{x, /?) = ±6'"'/^ e^2-^/^ (l){x, 0). (9) 

Since the partition function is left invariant by variable changes, we conclude that 

27TrT' 



Z[fii] = Z 



/"/ + 



r G Z. (10) 



This implies that ^ G [~;^) ;^]- 

At high temperatures, when the weak coupling approximation is consistent, a pertur- 
bative evaluation of the free energy F[fij] = — ln[Z]/f3 shows first-order phase transitions 
as a function of ^/ at fij = ^rw = (2r — 1)ttT/N (r G Z). These are the Roberge- Weiss 
transitions. At low temperatures, when the theories are typically strongly coupled, lattice 
studies suggest that F[fij] is a smooth function of m/T. The situation at very high and 
very low temperatures is therefore expected to be qualitatively as in Figure 1. The nature 
of the transition between these two behaviours depends on the detailed dynamics of the 
theories in question, as we have already discussed in the Introduction. 

2.1 Roberge- Weiss transitions at weak coupling 

At suitably high temperatures when the gauge theory is in a deconfined phase it is generically 
possible to compute the effective potential as a function of fij perturbatively [49]. The 
details of this calculation at the one-loop level are summarized in the Appendix. Here we 
state the result of this calculation. The basic idea involves computing an effective potential 
for the eigenvalues of the holonomy around the thermal circle (the Polyakov loop matrix), 
exp(i § A-r) in the presence of fundamental matter. The effective potential at one- loop is the 
sum of a gluonic (adjoint matter) piece and a flavour (fundamental matter) piece. While 



the gluonic contribution, Va, is ^Af invariant and independent of /i/, the terms arising from 
the flavour fields, Vf, are not Z^v invariant and explicitly depend on the imaginary chemical 
potential, 

K;flf ( {ai} ) = Va{ {aij} ) + Vf{{ai - m/T} ) , aij = ai - Oj. (11) 

Here e*"^, j = 1,2, ... N are the eigenvalues of the Polyakov loop matrix. Generically, at 
high temperatures, in the absence of fields in the fundamental of SU{N), Va will have N 
minima, each of which breaks the Z^r symmetry spontaneously. At these minima qj = 
2Trk/N for k = 0,1,2, . . . N — 1. Introduction of fundamental matter turns these into local 
minima with a single global minimum at a^ = for vanishing m. Since Vf depends only 
on the combination (oj — m/T), the global minimum will move to Qj = 2Tir/N when 
[iijT = 27rr/iV (r G Z). We conclude that the different local minima of the effective 
potential compete as /U/ is varied and there are first order phase transitions between two 
neighbouring minima whenever [iijT = (2r — 1)ti/N (r G Z). The perturbative result for 
the free energy across these transitions is computed in the Appendix. We quote the result 
here, specifically for a large-iV gauge theory coupled to both adjoint and fundamental 
matter fields. For such a theory with nj complex scalars and nj Weyl fermions, all in 
the fundamental representation of SU{N) and with mass niq, we find the free energy as a 
function of /x/ to be 



F[^iI\ = N■^^ 



rif /bosc ( ^ ) + nf /fcrmi ( ^ ) min ( -^ - — - ) , (12) 
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where 



and 



3 f^ 
/bose {x) = —^ y^ dy ^ z^- (13) 

2^ h sinh^ n^y-i + x^ 

3 /""^ o . 1 



/fcrmi {x) = -^ I y^ dy ^ — . (14) 

2^ ^ cosh^ (1^/2/2 + x2 



These functions are determined by the Bose-Einstein and Fermi-Dirac distributions for free 
particles. As nig/T — )■ 0, /bosc and /fcrmi approach 2 and 1 respectively. 

The free energy is a function with cusps (Figure 2) located at odd integer multiples of 
tt/N where phase transitions occur. It is a quadratic function only because we chose to look 
at the large- A^ limit which forces the Roberge- Weiss transition points to be closely spaced 
within intervals of 2tt/N. Hence only the second derivative of the free energy function at 
its minimum is relevant. The strength of the transition, as measured by the coefficient 
functions /bosc and /fcrmi varies with the ratio rUq/T. We have plotted these in Figure 2. 
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Figure 2: Left : The generic behaviour of the free energy across Roberge- Weiss transitions. Right: 
The functions /bose (blue) and /fcrmi (red) versus the ratio niq/T. 

Notable features are their behaviours for small and large masses, /bose is determined by 
the Bose-Einstein distribution which diverges for small masses and is responsible for the 
linear dependence on niq/T when niq <C T. For large masses (or low temperatures), both 
coefficients are exponentially Boltzmann-suppressed. 

We note here that at weak coupling only the ratio mq/T is relevant. Typically however, 
nonperturbative physics will generate another scale (such as Aqcd in QCD). In that case 
the notion of high and low temperatures and masses has an independent meaning and the 
weak coupling result can be trusted when T ^> Aqcd- In the model we study in this paper, 
the 't Hooft coupling A = Oym^ "^^^^ ^^ dialled to any fixed value (assuming the number 
of flavours is small compared to A^), and there is no dynamical scale. This means that the 
ratio rUq/T is the only parameter which can be varied at fixed 't Hooft coupling and we 
expect the above picture to be valid for A ^ 1. For large A, the perturbative evaluation is 
clearly not valid and we must use the dual gravity picture. 



3 AA = 4 theory coupled to A/^ = 2 flavours 

We will consider the A/" = 4 theory in flat space, and at finite temperature , coupled to A^ = 2 
supersymmetric matter transforming in the fundamental representation of the SU{N) gauge 
group. When the number, Nf, oi N = 2 hypermultiplets is small {Nf <^ N), the effect 
of these hypermultiplets on the dynamics of the A/" = 4 matter (such as the running of 
the gauge coupling) can be consistently neglected and the matter fields can be treated in a 
"probe" approximation. The AA = 4 degrees of freedom are always in the deconfined phase 
at finite temperature in flat space. This particular theory can be readily studied at strong 
coupling using the AdS/CFT correspondence as was first pointed out in [13]. 



3.1 The classical theory 

The AA = 4 theory can be naturally coupled to matter in the fundamental representation 
and thereby preserve a maximum of eight supercharges {J\f = 2 SUSY). In the language 
of AA = 1 superfields, the M = 1 hypermultiplet consists of two chiral superfields (Q, Q) 
transforming in the (iV, A^) representation of the SU{N) gauge group. The superfields in 
the hypermultiplet each consist of a squark and a quark {q,ip) and {q,tp). 

Introducing Nf hypermultiplets {Q^,Qi) and using the M = 1 SUSY notation, the 
superpotential of the theory with flavours can be written as 

W = ^ [Y.i^Q^'^sQ' + mgQi Q') + V2Tt{^3 [$i, $2]) I . (15) 

The three adjoint chiral multiplets of the AA = 4 theory are <I*i, <&2 and ^3. After coupling 
to the M = 2 hypermultiplets, <I>3 naturally becomes the scalar part of the M = 2 vector 
multiplet, while <5i^2 together constitute an AA = 2 adjoint hypermultiplet. Note that while 
it is possible to introduce distinct masses for each of the Nf flavours in question, we will 
only consider the case where all masses are equal. For equal quark masses the theory has a 
U{Nf) ~ U{1)b X SU{Nf) flavour symmetry group. The U{1)b is identified as the "baryon 
number" symmetry^ under which Q^ and Qi transform with equal and opposite charges. At 
the level of the Lagrangian, the mass term for the hypermultiplet induces the operator, 

C D — — (nig i^^i^i + mlqjq' + m^gjg* + \f2mqq\^:>,q' + +V2mqqi^zct'^ + h.c. j . (16) 
9ym ^ ' 

In what follows, the expectation value of the operator Sm = -q:^ will be referred to as the 
"quark condensate". 

3.2 D-brane picture 

At weak coupling (and zero temperature) in the D-brane picture in type IIB string theory, 
the M = 2 matter hypermultiplets can be introduced by placing N j coincident D7-branes 
in the presence of A^ coincident D3-branes. The low energy dynamics of the latter is AA = 4 
SUSY gauge theory. The D3 and D7-brane worldvolumes span the coordinates indicated 
below in the ambient ten dimensional spacetime, 

Coordinates xq x\ X2 X3 X4 X5 xq xy xs xg 
D3 X X X X 

D7 xxxxxxxx 



^For massless quarks the theory has the global symmetry group, U{1)b x SU{Nf) x SU{2)r x U{1)r x 
5'C/(2)$. The SU{2)ji x U{1)r factor is the R-symmetry group of the N = 2 gauge theory. 



The fluctuations of the open strings stretching between the D3 and D7-branes yield the 
Nf hyperniultiplets transforming in the fundamental representation of the SU{N) gauge 
symmetry of the D3-branes. These also transform under the U{Nf) ~ U{1)b x SU{Nf) 
symmetry group associated to the D7's, which appears as a flavour symmetry at low energies. 
The separation in the xs-xg plane, between the stack of D3-branes and the flavour branes 
translates to a mass niq for the hyperniultiplets. For a separation i in this plane, the 
hypermultiplet mass is given by the string tension times the separation, niq = 2:;^- 

4 Thermodynamics at strong coupling 

In the large- iV and strong 't Hooft coupling limit (A = Qym^ ^ 1)' ^^^ fundamental 
hyperniultiplets are introduced as probe D7-branes into the geometry sourced by the N 
D3-branes, namely AdS^ x S^ [13]. At finite temperature, in the deconfined phase, the 
corresponding geometry is the Schwarzschild black hole in AdS^ x S^ (with fiat horizon). 
Wick rotating to Euclidean signature, the spacetime metric is 

ds^ = {j{r)dT^ + ;^ + i^^' + ^' ^^s) , (17) 

where 

f{r) = ;^ (^' - ^) , rH = ttTR^ R^ = M9sN)a'^ = g^^N a'" , (18) 

and R is the AdS radius. 

In the Euclidean picture the geometry ends smoothly at r = r//, provided we make 
the periodic thermal identification t c:::^ t + 1/T. Topologically, the Euclidean black hole 
spacetime is M.^ x D2 x S^ . 

As suggested by the weak coupling setup in Section (3.2), the flavour D7-branes must 
wrap an S*^ C S^ while filling the AdS^ directions [13,15]. The embedding of the D7-brane is 
characterized by the "slipping angle" which determines the size of the three-sphere wrapped 
by the D7-brane. The slipping angle 6 is defined by the following paramerization of the S^ 
metric 

dnl = de^ + sm^ed(l)^ + cos'^edQl; 0<9<^. (19) 

It is more convenient to rewrite the background metric in a FefFer man- Graham type coor- 
dinate system [50,51]. We use the conventions of [21], to write the black hole metric in the 
form 

d'' = lT^{jdr' + h dx'^ + ^ {dp' + p' dnl) , (20) 

10 



4 4 

h{p) ^1-^, h{p) ^1+^. 

P P 

This form of the metric is obtained after the coordinate change r^ = ^ (p^ + rfj/p^), and 
we take the metric of the five-sphere as in Eq.(19). The shpping angle 6 will be dual to the 
quark mass operator S^, and depend only on the radial coordinate of AdS space, 

= 0{p). (21) 

This parametrization is convenient for making contact with the weak coupling D-brane 
picture. Focussing on the six directions transverse to the D3-branes, we write 

dp^ + p'^dnl = dpi + pI d(p'^ + dpi + pI dQJ , (22) 

with 

pi=pcos6, p2 = psm9. (23) 

The xs — xg plane can now be identified with the {p2, (f>) plane. We choose the D7-brane to 
he at (j) = 0, i.e. the xg axis. The asymptotic value of the D3-D7 separation, given by p2 
for large p, yields the hypermultiplet mass niq. 

4.1 Zjv breaking in thermal N' = 4 theory 

As argued in [52], one may understand the spontaneous breaking of the Zjv center symmetry 
at high temperatures as follows . We first note that the disk D2 in the black hole geometry 
naturally allows turning on a Neveu-Schwarz potential B with dB = and 

This yields a phase e*" for the Polyakov-Maldacena loop^ [53] computed by wrapping a 
fundamental string worldsheet on the cigar D2. Thus a is defined modulo 27r. One then 
finds that an effective potential for a is induced if we consider terms in the type IIB effective 
action that depend on the RR 3- form i^(3), 

^IIB ^ (^ (/ d''x ^g ^F(|) - j F(5) A F(3) A s) . (25) 



^Strictly speaking, the Wilson loop computed by the wrapped string world-sheet is the supersymmetric 
version of the standard Polyakov loop and includes the adjoint scalar fields of A/" = 4 SYM in its definition, 
U = expi ^(At— i^^(r)$^). This is an important distinction in principle, but will not influence our discussion 
below, as the Polyakov-Maldacena loop has the same transformation properties under Zjv transformations 
as the standard Polyakov loop order parameter. 

11 



Integrating over S^ x D2, we may obtain the effective theory for the ^(3) field on W^. Since 
F(3) has only one non-zero component -F123 on M^ we have 



Ses[Fi3)] = J d\ 



(^123) - 73^7^^123 



T^2^^^lV ''■" 47r2/r 



(26) 



where we have used the fact that the type IIB background has N units of five-form flux 
through S*^, /cs F^/{2t:Is)* = N . Adapting the arguments arising in the context of U{1) 
gauge theory in two dimensions to the present situation [52] , gauge invariance of quantum 
states implies the quantization of -Fi23- States are characterized by a constant (quantized) 
value for the field strength (similarly to constant electric fields in two dimensions) and one 
finally finds the effective potential for the phase of Polyakov-Maldacena loop to be 

1 / o„ \ 2 

VA{a) = Att'^N^T^ - min ( a - — r I . (27) 

A rez y ^ J 

The potential is minimized when a = liir/N for r € Z and this leads to a spontaneous 
breaking of the Z^v symmetry of the N = A theory in the deconfined, high temperature 
phase. Note that this effective potential is 0{N'^)^ as expected for a leading order large- A'^ 
effect and proportional to T^, as expected from dimensional analysis. Interestingly, it is 
also proportional to 1/A in the strongly coupled theory. 

4.2 The effect of flavours 

Now we would like to calculate how the effective potential for a is modified upon introducing 
fundamental flavours and in the presence of an imaginary chemical potential for quark 
number. The basic effect of introducing flavour D7-branes was already addressed in [54] 
and was argued to lift the degeneracy of the N vacua discussed above. The overall dynamics 
of Nf flavour D7-branes, ignoring the non-Abelian world-volume degrees of freedom in the 
probe limit {Nf <^ N), can be described by the (Euclidean) Dirac-Born-Infeld action 

Sd7 = Nf Td7 I d'^ie-^ Vdet(*(/ + (27ra'F + *B)) , (28) 

where *g and *B are the pullbacks of the spacetime metric and the Neveu-Schwarz two- form 
potential on to the D7-brane worldvolume. In addition we have that the dilaton e~* = 1 
in the background dual to the AA = 4 SUSY gauge theory. The D7-brane tension can be 
written in terms of gauge theory parameters A, A^ (and R, the AdS radius) as 

^^^ ~ gsa'S2ny ~ SItt^R^' ^^^^ 
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We will take the D7-brane probe to wrap the S^ in Eq.(19) at an angle 6 and fill 
the directions r, x and p. For the most part, we will primarily be interested in classical 
configurations that completely wrap D2, the cigar in the Euclidean black hole, corresponding 
to "melted meson" states [18,21]. For these configurations it will be consistent to turn on 
a field strength F along the world-volume. A real value for F in the Euclidean action 
will correspond to an imaginary gauge potential upon Wick rotating back to Lorentzian 
signature. 

4.3 Imaginary chemical potential 

In the boundary field theory formulated on Minkowski space, a chemical potential for a 
global f/(l) symmetry can be thought of as follows: We imagine that the global symmetry 
is gauged and then turn on a constant background value for the time component of the 
gauge field, while setting other components to zero. A real chemical potential is obtained 
by assigning an imaginary constant value for this background gauge field in the Euclidean 
formulation. Equivalently, to get imaginary chemical potential for the U{1)b global sym- 
metry, in the Euclidean theory we must introduce a real value for the time component of 
the background gauge field associated to U{1)b- This is very natural in the AdS/CFT 
framework. 

The gauge field associated to the overall U{1) factor describing the center-of-mass dy- 
namics of the probe D7-branes, is dual to U{1)b- In the Euclidean theory, the imaginary 
chemical potential m corresponds to a constant real boundary value for the gauge field Ar 
on the world-volume of the D7-brane, 

p,i = -T dr dpFrp= lim Arip). (30) 

Jo Jth ^^°^ 

Here we choose the gauge freedom to set Ap = {) and we require that ^^-(r//) = in order 
to avoid a singularity at the tip of the cigar. We are implicitly assuming that the D7-brane 
wraps the black hole cigar D2 completely. The case where it turns back before getting to 
the horizon will be discussed separately. 

The action for the D7-brane probe configuration described above is 



Sm = 2TT^NfTD7 [pd^xTdpil^hh cos^ oJl + p^e'{p)^ + 2(27ra'F,p + *Brp?^ 

(31) 
We note that in the absence of a S field in the background, this is the same action used 
in [21], provided we make the replacement F — )■ iF. The qualitative nature of the family of 
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solutions is, however, different to that encountered in [21] for real baryon number chemical 
potential. 

4.3.1 Periodicity in /i/ 

The periodicity in fij discussed earlier is naturally built into the supergravity dual because 
the D-brane action only depends on the combination B + Ina'F. This means that the 
classical solutions to the DBI equations of motion will be charactarized by the combination 

A shift of /i/ by 2itT can be absorbed by a — t- a + 2tt. Hence m is only defined modulo 
27rT as expected from general considerations. 

As at weak coupling, there is a further periodicity visible in the finite temperature theory 
at strong coupling. The shift /[// — )■ ^/ + 2itT/N can be compensated by a corresponding 
shift of a by the same amount, leaving the physics unchanged since the "glue" sector of 
the theory (M = 4 matter) is invariant under such shifts of a (from Eq.(27)). Note that 
this "periodicity" involves a jump from one Z^v vacuum to another. Therefore, this system 
should be expected to display the infinite sequence of Roberge- Weiss transitions as the 
effective potential will be a sum of two pieces, one of which is Zjv invariant, while the other 
piece, induced by the fundamental flavours, depends only on the combination a — fJ,i/T. 

4.4 Solving the DBI equations of motion 

The solutions of the D3/D7 system at finite temperature, with and without a real chemical 
potential, have been extensively studied in earlier works [15-23,25,31]. Our aim, however, 
is to study the solutions arising in the presence of the imaginary chemical potential fij ^ 0. 
These will be distinct from the ones found, for example, in [21,22,25]. Specifically, for a 
real chemical potential, the term under the square root in Eq.(31) would be of the form 



1 + p'^9'^ — 2{2Tra' F)^h/h'^, yielding solutions with different behaviour to the ones we will 
find below. 

From the solutions to the DBI equations of motion following from Eq.(31), we will 
reconstruct the dependence of the effective potential on a — fij/T, for generic quark masses. 

We work in the gauge Ap = 0, so that F^p = —dpAr- Since the action (31) does not 
depend explicitly on A^, the conjugate momentum is conserved. In the usual situation 
with real chemical potentials, the variable conjugate to A^ is the quark density. We will 
continue to refer to this conjugate variable as the "charge density" d, in the situation with 
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an imaginary chemical potential, 

— a'lr = ^^ ^/^^7(27rQ; ) — — cos^fel ^ '^ '^ — _ . (33) 



^Jl + p^9'^ + 2{27ra'Frp + Brpy-^^ 
This is a constant of motion. It is useful to define the dimensionless combination 

d ^ d(2W)-i {27:'NjTo7rj,)-' = ^7^^^- (34) 

The solutions can be characterized by this dimensionless density which, as we will see 
below, is restricted to take values less than unity. Solving algebraically for the combination 
{2TTa'F + B), using the above equations we obtain 



B.P + 2WF., = d ^^^^^nrww _ (35) 



h(h^Ccos^e-sd? 



3-' 

' H 



For solutions that reach the horizon of the black hole (i.e. which wrap the cigar D2), we 
deduce from the above expression that 

|J| < cos^^l^^,^. (36) 

The explicit relation between the density d and the chemical potential /U/ follows from 
integrating Eq.(35) 



''^ W (1 + y-4) (^(1 + y-4)3y6 cos6 e - Sd?\ 

Here, the dimensionless variable y = pjru. 

The D7-brane configurations that reach the horizon are referred to as "black hole embed- 
dings" . Analysis of the spectral functions of various fluctuations around such solutions re- 
veals a continuous, gapless spectrum [55,56] leading to the interpretation that these describe 
a phase where the "mesonic" fluctuations have melted in the high temperature plasma. 

4.4.1 The action 

Evaluated on the solution for B + 2'koIF ^ the unrenormalized D7-brane action is 



^(l + y-4)3y6cos6^_8j2 



Note that this action cannot be varied to obtain the correct equations of motion for 0{y). 
The equation of motion for 9 follows as usual by varying (31) or by using a different action 
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where we trade Aj- for F^-p as the independent variable, the equation for the latter being 
algebraic. This is achieved by the Legendre transform, 

Sd7 ^ SJ^T = Sdt - d j P d^x I dp (f^p + ^\ . (39) 

The resulting Legendre transformed action can be evaluated on the solution for F, to yield 

7^3 [■ roD / Rrl'^ \ ^^"^ 

5B? = AiV/iV— / d^x / dy y' hh y^TTWW 1 " ITI^ r- cos^ 6 . (40) 

b^ J Ji \ h-^ y° cos° t^/ 



It is easily verified that the corresponding expressions for real densities and real chemical 
potential can be obtained by the replacement d — )■ id. 

4.4.2 ruq = 0: constant solutions 

It is straightforward to see that the DBI action admits constant solutions with ^ = for a 
range of densities d. Since the constant solutions extend all the way to the horizon, they 
describe a phase of melted mesons. They also correspond to having vanishing masses for 
the fundamental hypermultiplets. In general, the large p asymptotics of 6(p) can be easily 
seen to be 

<'W = ^ + % + .... (41) 

P p'^ 

where the coefficients 6m\ and 9^2) ^^^ related to the quark mass and condensate respectively. 
According to the AdS/CFT dictionary, this boundary behaviour signals the fact that 9 is 
dual to a dimenson 3 operator, namely, the quark bilinear. From the earlier description of 
the D3-D7 system and from Eq.(22), the asymptotic value of p sin 9 is the D3-D7 separation, 
which is the quark mass i.e., 0(o) = 2TTa' niq. The coefficient of the subleading term, 
namely 0(2), determines the VEV of the same operator, i.e the quark condensate (Sm.) = 
(^i^i) + •••• 

It follows then that the solution with 9 = has niq = 0, (S^) = 0. When the slipping 
angle 9 is vanishing, the D7-brane wraps the equatorial S^ inside the five-sphere in the 
geometry. For this solution, the contribution of the flavours to the effective potential for a 
can be computed analytically^. From Eq.(37) we find 



1 f pi\ d f°^ ^^ 1 ^,J,-2/3o /^l 1 

\/y^ - d? 



Txy-T)--2h ^^:7==-T^'^'"^^^^H3'2)' (^^) 



*The same expressions were obtained in [25] for the zero mass 'black hole embedding' solutions for real 
baryonic chemical potential. 



16 



where Bz{a,b) is the incomplete beta function and the integration variable y is related to 
y as y"^ = 2 (y^ + 2/"^)- The action is formally divergent and requires renormalizaton which 
is achieved by a simple subtraction for the constant solutions, 



cren 
'-'D7 



XNjN^J,'.^^ 



dy 



y 



VW^dP 



y 



(43) 



rp3 

XNfN — 



'f 



96 



d^x 



\d\'/'B,. 



d? 



2 1\ 3 
'3'2 j ^2 



The effective potential induced by the flavours can be found by eliminating d from equations 
(42) and (43). It is important to note that since a is a phase defined modulo 27r, and likewise, 
the imaginary chemical potential ^u/ is also a periodic variable, the combination a — ^i/T 
is bounded. 




-0.1 0.1 



Figure 3: The potential for a from the flavours, as a function of [a — /i//r)/vA. The curve in 
green is the full potential implied by the DBI action. The curve in blue is the quadratic piece that 
is consistent to keep in the A -^ 00 limit. The normalization constant M — NfNT'^X/16. 



In the strong coupling limit A ^ 1, we must therefore have that 

J=-^(a-/i,/r) + 0(I/A), 



(44) 



by inverting Eq.(42). Since |(i| ^ 1 at strong coupling, we can expand the action (43) for 
small d to obtain the contribution to the effective potential from the flavours, 

2 



Vj = NfNT'[a-^) +0{l/V\). a,^G[-vr,7r], 



(45) 



which is independent of A at the leading order in the strong coupling limit. The full effective 
potential at strong 't Hooft coupling with Nf <^ N in the deconfined phase is then 



^2 ]\t2 



rp4 



V,s =VA + Vf = mm Att'N'— [» - ^ ) + ^jNT'' [a ^ ^ 



27rr\ ■ 



^v 



(46) 
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From this effective potential it is obvious that the degeneracy of the Z^v vacua has been 
Hfted by the massless fundamental flavour fields. Curiously, while the flavour contribution 
is 0{\/N) suppressed compared to the leading term in the large A^ limit, it is unsuppressed 
by powers of the 't Hooft coupling A in the strong coupling limit. The order of limits is, of 
course, unambiguous in the present context since A^ is taken to infinity first, keeping the 't 
Hooft coupling fixed and large. 






Figure 4: The full effective potential as a function of a from the flavours for three different values 
of fJ-i/T. As fJ-i/T is increased from left to right, the local minima of the effective potential compete, 
resulting in a first order phase transition. 

There are N local minima, but a unique global minimum at a = when ^/ = 0. As 
/// is increased smoothly from zero and approaches fij = ir/N, two neighboring minima 
become degenerate and any further increase in /// results in a first order phase transition 
from the a = vacuum to the vacuum with a = 2tt/N (see Fig. (4)). 



Free energy: As the phase transitions at fii/T = (2r — 1)-^ for r G Z are first order, 
it is useful to compute the free energy as a function of the imaginary chemical potential ///. 
The free energy (the grand potential) density can be calculated by minimizing the effective 
potential Eq.(46) with respect to a for all /i/ and we find 



f[n 



A>1 



NfNT'^ ( 1 + 



A A^ 



~ NfNT'^ min ( ^ - ^^ 



47r2 A^ 
fij 2'Kr 



mm 



T 



2'Kr 



(47) 



r&Z \T N 

While the quadratic behaviour is simply a consequence of large A^, it is interesting to 
compare the coefficient with the result at weak coupling and zero hypermultiplet mass 
Eq.(12). Since we have 2Nf complex scalars and 2Nf Weyl fermions, the free energy at 
weak coupling (and ruq = 0) is 

2 



f[n 



A<1 



1 MAT ^4. ■ I f^i 27rr 
- JN JS fl mm 

2 ^ r&\T N 



(48) 
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We observe that the strong and weak couphng resuhs differ by a factor of 2. 

4.4.3 ruq / non-constant solutions 

Solutions with non-zero quark mass also happen to be non-constant and these need to be 
understood numerically. For smooth solutions that get to the horizon from the boundary, 
we require that the S^ C S^ wrapped by the D7-brane remains non-vanishing along the 
solution. As is well-known, there are two classes of non-constant D7-brane solutions: those 
that get all the way to the horizon (the so-called "black- hole embeddings" ) , and those that 
end smoothly before they get to the horizon. For both classes of solutions, the quark mass 
and the condensate can be read off from the asymptotic behaviour (41) near the boundary of 
AdS^. Specifically, in terms of the dimensionless variable y = p/rn, the boundary behaviour 
of the slipping angle is 

Note that it is only the ratio niq/T which is relevant here, since the theory without hyper- 
multiplets is the conformal M = 4 theory without an intrinsic scale. We thus define the 
natural dimensionless ratio that characterizes these solutions, 

T Vat' ^ ' 

where the mass scale M = 2mq/y\ is the characteristic scale of meson bound states at zero 
temperature (and strong coupling) [30]. 

The Euclidean action for these solutions requires careful renormalization via subtrac- 
tions to yield finite results. The correct procedure of holographic renormalization has been 
developed for probe Dp-branes in AdS spacetime [31,51]. The result of this procedure in 
the present context is that the renormalized action is defined as an integral over y from 
y = 1 (the IR) to y = A (a UV cutoff) along with counterterms defined locally at y = A, 

S'3^ = SdiW - AiV/iV^ Qa^ - \K'e{Kf + ^A'^(A)') • (51) 

The first two counterterms are divergent while the third is a finite subtraction. Further 
counterterms are necessary if the boundary is not fiat. 

Black-Hole embeddings: Solutions extending all the way to the black hole horizon, of- 
ten referred to as "black-hole embeddings" of the D7-brane, satisfy 9'{y = 1) = and 
< 0{y = 1) < ^. Thus for these embeddings, the S'^ remains blown-up all along the 
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solution, while the S^ (the thermal circle) shrinks at the horizon, so the D7-brane caps off 
smoothly. By varying 6{y = 1) we can explore different values of the hypermultiplet mass 
M/T. These Euclidean solutions can exist only if 



\d\ <cos3 0|j^. 



I- 



(52) 



The black hole embeddings describe a high temperature phase where mesons have melted. 
For a fixed density d, solutions with different hypermultiplet masses can be found by 
dialling 9{y = 1) provided the condition (52) is satisfied. As in the massless case analyzed 
above, we will obtain (numerically) the effective potential as a function of (a — /i//T)/vA. 
We summarize the properties of the solutions below. 

• The solutions for finite density < \d\ < 1, are all of the type described above and are 
qualitatively similar to black-hole embeddings found earlier for vanishing baryonic 
chemical potential [15, 18]. They are qualitatively distinct from the configurations 
with real chemical potential [21,22,25,26]. We elaborate on this below. 

• For any given density \d\ (between and 1) there is a fixed value of M/T, above 
which solutions cease to exist. The allowed range of masses decreases with increasing 
\d\ (Figure 5). When the density d approaches unity, all allowed solutions bunch up 




Figure 5: The values of ^ and the density variable d (conjugate to imaginary fi) for which the 
black hole embeddings or melted meson solutions exist. 

near the zero mass constant solution. This is illustrated in Figure 6 where we plot 
solutions in the pi — p2 plane (p2 = psin^ and pi = pcosO). In contrast to this 
picture, for a real chemical potential, it was observed (e.g. in [21, 25]) that black 
hole embeddings can exist at all possible values of the ratio M/T. For large values of 
the mass (or low temperatures), these solutions could be viewed as a probe D7-brane 
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with strings (quarks) attached which drop into the black hole horizon. The absence of 
these kinds of solutions substantially alters the phase diagram for imaginary chemical 
potential. 
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Figure 6: D7-brane solutions extending to the horizon for d = 0.01, 0.5 and 0.98. For each density 
we have plotted a family of solutions with increasing masses up to the maximum value allowed for 
that density. The asymptotic distance from the pi axis gives the ratio M/T. The black hole horizon 
is the blue curve pi + P2 = ^h- 



2m„ 



The largest value of the ratio ^ = ^y=^ for which the black hole embedding exists is 
w 1.30 and this occurs when d = and is the critical embedding solution (Figure 10). 
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Figure 7: The action as a function of M/T for d = and d — 0.5. 

• By evaluating the action and the expectation value of the quark condensate operator 
(Sm), for a given density d, we may see that for large enough masses, there exist 
two classical solutions with the same mass but different values for the action and 
condensate (Figure 7). This feature, well-known at zero density, also appears for non- 
zero densities conjugate to an imaginary chemical potential. The solution with the 
lower action is then chosen to be the physical one. 

The effective potential: The contribution to the effective potential for the phase of 
the Polyakov-Maldacena loop a, can be found numerically (Figure 8). However, we can 
already be fairly precise about the form of this contribution due to Eq. (37): since a and 
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Figure 8: The effective potential from flavours for increasing values of mq, for M/T = 0, 0.24, 0.48, 
0.72, 0.96 and 1.2. The potential tends to flatten out as M/T approaches the critical value of 1.3. 

/X//T are bounded, we must have that \d\ ~ ©(l/vA) ^ 1. Notice that this scaling of 
d, actually implies that the dimensionful density d ~ NfNT'^, which is a natural scaling 
for Nf quark flavors. In taking d ~ ©(l/vA), we are also assuming that the integral in 
Eq. (37) is well-behaved as d — )• and we can check numerically that this is the case 
for all hypermultiplet masses in the range < ^ < 1.30 corresponding to the black hole 
embedding solutions. In addition since the D7-brane action is a function of d^, for small d, 
it can be expanded in powers of d^. In the large A limit, the quadratic piece alone dominates 
and therefore we find that the effective potential is of the form 

V.^=Va + V, = min4^2iV2^(a-^)' + iV,iVTV(f ) («-f )'• (53) 



2r?i„ 



Here M is the mass scale —t^- and / is a dimensionless function of -jr- The coefficient 
/ ( jt) is unity when mq = and decreases monotonically with increasing mass, eventually 
approaching ~ 0.2 as -^ is dialled to its maximum value of approximately 1.3 (Figure 9). 
We should emphasize that in deriving the effective potential above, in the large A and large 
N limit, we only needed to use solutions with small d, strictly only those in the vicinity of 
d ~ 0{l/v\). When we turn to the low temperature regime below, we will see that this 
is actually crucial for ensuring a consistent picture of the transition between high and low 
temperatures. 

To summarize, for all values of the parameter -^ for which black hole embeddings of the 
D7-brane exist, with an imaginary chemical potential, first order Roberge- Weiss transitions 
will occur at ^ = (2r — 1);^, r G Z. This means that, for fixed hypermultiplet mass m-g, 
the RW transition will occur for temperatures T > 2mq/(1.3vA). We will determine the 
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Figure 9: The dimensionless coefficient / in the quadratic effective potential Vf due to hypermul- 
tipltes. 

value of this temperature, representing the RW endpoint, more precisely below. 

The monotonic decrease in f{M/T) with decreasing temperature (or increasing mass 
ruq) is in accord with intuition from weak coupling. The quadratic potential for a can be 
interpreted as a thermal contribution to the mass of the mode a from the flavours. That this 
thermal correction should decrease as the flavours are made heavier (eventually decoupling 
for infinite mass), appears intuitively to be consistent. At weak coupling, it would be natural 
to identify the mass of the mode a with the Debye mass (inverse electric screening length). 
It is unclear whether any such interpretation should be possible at strong coupling. Finally, 
we can again write the free energy as a function of ///, (taking XNf/N ^1), 

F[^,J] = NfNT' f (^) mm (a - f,j/Tf. (54) 



mi 



Given our numerical results for /, we can say that near niq = 0, / ~ 1 — f"{^) o\^i ■ In 
contrast, in the weakly coupled theory the coefficient in Eq. (12) (with Uf = hf = 2Nf) 
decreases linearly with mass, for small enough mass, and is independent of the weak 't Hooft 
coupling at leading order 

1 1 777 TT 

g ( /bose(m,/T) + /fermi(m,/r) ) ^ - - LA _ + , , , (55) 

The origin of the linear term can be traced to the pole in the Bose-Einstein distribution for 
free massless bosons at zero momentum. The absence of a linear term at strong coupling 
is perhaps indicative of the strong dressing undergone by the perturbative scalar degrees of 
freedom. It is also worth noting that at strong coupling, the natural dimensionful parameter 
is the meson mass scale M, rather than the hypermultiplet quark mass rriq. 

Low temperatures/Large masses, Af > 1.3: The critical black hole embedding of 
the D7-brane has the property that the slipping angle 6 of the solution approaches ^ at the 
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Figure 10: The well known transition at zero density, from black hole embeddings (falling into the 
black hole) to Minkowski embedding probe D7-branes placed away from the horizon. 

horizon and the corresponding 5^ wrapped by the probe brane shrinks. Hence, the only 
solutions with masses higher than ^ ~ 1.3 are those where the S^ wrapped by the D7- 
brane shrinks before the D7-brane gets to the horizon. For these embeddings the thermal S^ 
remains non-vanishing while the S^ shrinks smoothly before the brane gets to the horizon. 
Since the S'^ shrinks at the tip of such a solution, we are forced to have d = due to Eq.(52) 
and hence 

Frp + 2Tra'Brp = 0. (56) 

The D7-brane action has no dependence on either a or ///. Hence there is no potential 
induced by flavours for the phase a of the Polyakov-Maldacena loop: Vf{a — fij/T) vanishes 
identically and the physics is completely smooth as a function fij. Note that despite the 
fact that the field strength F is fixed by B, we can always add a constant to the world- 
volume gauge potential Ar, which will not change the field strength, and can be interpreted 
as a chemical potential fij in the boundary field theory. Hence, in the grand canonical 
ensemble, at low temperatures (T <^ M), for every /i/ the only solution possible is the 
d = Minkowski embedding representing unmelted mesons. The situation is quite distinct 
from the case of real chemical potential [26] wherein, at fixed low temperatures, there is a 
transition from Minkowski embedding solutions to the spiky black hole embeddings as the 
chemical potential is increased. 

The solutions with d = are the same as those originally obtained in [15, 18,20]. The 
transition from the black hole embeddings to the second class of solutions involves a topology 
change: the first category of solutions exhibits a shrinking S^ and a finite S"^, while the 
second has the S^ shrinking in the interior. The transition between these is well-known and 
is a first order phase transition as may be inferred from plotting the action as a function of 
M/T (Figure 11) for the solution with d = 0. The phase transition between the two classes 
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Figure 11: Left : The action for d — 0, D7-brane embeddings. As the ratio M/T is increased, just 
before the black hole embedding (melted meson) ceases to exist, there is a transition to the family of 
so-called "Minkowski embeddings" corresponding to unmelted mesons. Right: In the second figure 
we see that the action (grand potential) for black hole embeddings with d > 0, (0.004 (orange) and 
0.008 (blue)), is always larger than the action for the d = solutions (black). There is no Minkowski 
embedding solution with d ^ 0. 

of solutions is a "meson-melting" transition. At low enough temperatures (or large enough 
hypermultiplet mass), the only allowed D7-brane embeddings are the unmelted mesons with 
zero density d = 0. The spectrum of fluctuations about this solution exhibits a mass gap 
and a discrete spectrum [30] at zero temperature. The mass gap in the meson spectrum is 

771 

Mgap = 47r^ = 2ttM. (57) 

vA 

We also see that as the density d is increased from zero, for small enough densities 
d <^ 1, the black hole embeddings with d ^ coexist with the zero density d = Minkowski 
embedding solution for a range of temperatures (Figure 11). This will be important for 
determining the phase diagram and the shape of the phase boundaries in the grand canonical 
ensemble. With increasing density, there is a critical value of d, (> 10^^) beyond which the 
curves for black hole embeddings (as in Figure 11) cease to intersect the Minkowski solution 
at d = 0. They turn back and terminate at progressively smaller values of M/T (Figure 7). 
The presence of such configurations in the grand canonical ensemble could make the low 
temperature region of the phase diagram inaccessible. However, it should be clear from the 
preceding discussions that these configurations are essentially removed as a consequence of 
the Roberge- Weiss transition. The RW transitions occur at values of {a — fJ-i/T) that are 
parametrically suppressed by 1/A^ corresponding to parametrically small values of d. The 
configurations with such low densities will always be of the kind depicted in Figure 11. 
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5 Phase diagram 

The above analysis of the dominant configurations at different densities and temperatures 
now ahows us to determine the phase diagram of the theory as a function of the temper- 
ature/meson mass scale ratio T/M and the imaginary chemical potential /x/. We have 
already established that at low temperature the physics is independent of /i/ and q, the 
phase of the Polyakov loop, and is dominated by the unmelted meson (Minkowski em- 
beddings) configurations with vanishing density d = 0. In addition, we have seen that at 
high temperature the dominant configurations are the black hole embeddings representing 
melted mesons with d 7^ 0. For these configurations the system experiences a first order 
phase transition at fii/T = (2r — \)tt/N (r € Z), characterized by a discrete jump in the 
phase of the Polyakov loop. Specifially, a jumps from 27r(r — 1)/N to 2Tir/N. 

What remains is to understand the phase boundary between the melted and unmelted 
meson phases as a function of /i/. We first recall that in the melted phase d is proportional 
to {fJ-i/T — a)/vA, for small enough d. As depicted in Figure 11, in the grand canonical 
ensemble, the melted meson configuration with small d ^ will then have to compete with 
the d = Minkowski embedding as the temperature is decreased. This results in a first 
order transition between the two phases. It is clear from Figure 11 that as the density 
increases (for positive d, say), the transition temperature increases as well (equivalently 
M/T decreases). The dependence of this transition temperature on the chemical potential 
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Figure 12: The melting transition temperature at different values of /i/ and fixed a. 

determines the shape of the phase boundary in the /U/ — T plane. We can numerically 
obtain the melting transition temperatures Tc{fj,i) for different values of /// (corresponding 
to different densities d in the melted phase). The result is shown in Figure 12. For 

«=^, (2r-l)-< ^ < (2r + l)-, r G Z, (58) 
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we find that 

T.ifij) - To = K ^ (-^^ - a] +..., K^33.5 (59) 

where Tq = Tc(0) ~ 0.77M. We have only kept terms to quadratic order on the right hand 
side for two reasons: First, since the difference /U//T — a is bounded, higher order terms are 
suppressed at strong coupling (A ^ 1). Secondly, we also have that {fJ-i/T — a\ < tt/N, so 
higher order terms are suppressed by powers of 1/N'^. At leading order in 1/A, Tc{jxi) may 
therefore be replaced by Tq on the right-hand side of (59). 

The first-order phase boundary between the melted phase at high temperature and the 
unmelted phase at low temperature has the curvature as shown in Figure 1. Combining this 
result with the Roberge- Weiss lines at /U/ = /Urw yields an infinite set of points at which 
three first order transition lines meet. At each of these triple points, three phases labelled 
by distinct values of the Polyakov loop coexist. The location of the triple points in the 
m — T plane is given by (^rw^Trw); where 

^ = ('^-')iV' ^ = 1+AA^' ^^°) 

with r G Z and n = KM/T ~ 43.5. The iV-scaling of the formulae may suggest that these 
should be subleading effects in the large- A^ limit and therefore not consistently incorporated 
in a classical supergravity approximation. However the factors of 1/A^ arise from purely 
kincmatical considerations, namely the symmetry of the theory under shifts ^i — )■ ij,i+2it/N. 
When A^ is large, the RW transitions are closely spaced, resulting in the A^-dependence 
in the Roberge- Weiss temperature. The curvature of the phase boundary in Figure 1 is 
independent of A^, see Eq. (2). 

5.1 Real chemical potential and analytical continuation 

For real values of the quark density a rich family of solutions to the DBI equations of motion 
exists, after making the replacement F ^^ iF (and setting i? = 0) in Eq. (31). In particular, 
black hole embeddings, representing melted mesons, exist for all values of the ratio M/T 
in the canonical ensemble with fixed quark number density. These include low temperature 
(M ^ T) solutions which can be viewed as Minkowski embeddings attached to a spike 
consisting of strings (quarks) falling into the black hole horizon. 

The phase diagram for real chemical potential in the grand canonical ensemble [22,24-26] 
bears little relation to Figure 1, either at low or high temperatures. With a real quark 
chemical potential, there is a single line of phase transitions separating a zero density phase 
from the dissociated or melted meson phase. For small chemical potential the transition 
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line is first order, while at larger chemical potential the transition is expected to be a 
continuous one [25,57,58], terminating at /i = ruq. In fact it has been argued in [57] that 
for low temperatures the transition line is actually third order and connects up with the 
first order line at a tricritical point. 

In the context of QCD, the main motivation behind the exploration of the phase diagram 
as a function of imaginary chemical potential has been to determine the phase structure for 
real chemical potential, via the use of analytic continuation from imaginary ^ = ifj,j (with 
fj? < 0) to real n (with fi^ > 0) [34-43]. It is therefore natural to ask what aspects of the 
phase diagram can be captured by analytic continuation in the holographic model we studied 
here. This is particularly interesting given the qualitative difference in the nature of black 
hole embeddings and the D7-brane solutions at fixed density. At large real or imaginary 
chemical potential the phase structures are manifestly different due to the presence of the 
Roberge- Weiss transitions. However it is expected that the curvature of the meson melting 
line for imaginary chemical potential, with /i^ < 0, is directly related to the curvature of the 
first order transition line for real chemical potential, with fi^ > 0. Repeating the analysis 
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Figure 13: The line of first order meson-nielting transitions for real chemical potential, /i > 0. 

for real chemical potential we find that the curvature of the first order line for /u^ > is 
given by (see Figure 13) 



T,(/.)-ro«-33.5y^, (61) 



^0 

matching up with our result for imaginary [i. It is reassuring that despite the differences 
between individual D7-brane probe solutions for /i^ ^ 0, the theory is analytic in /u^ near 
^2 = 0. 
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6 Discussion 

In the presence of an imaginary quark chemical potential fij the D3-D7 holographic setup 
at strong coupling displays a first-order meson melting transition at T = Tc{^i) and a set 
of first-order Roberge- Weiss transitions at fixed /U/ = ^yon in the high-temperature phase 
T > Tc{^i). The first-order lines join at a series of triple points in the fij — T plane. We 
have determined the location of the triple points and the curvature of the phase boundaries, 
and confirmed that the theory is analytic in ^u^ near fi^ = 0. 

Since the adjoint degrees of freedom are always deconfined in this model, the quark (or 
meson) mass and the temperature cannot be varied independently in the absence of another 
scale. The model of Sakai and Sugimoto [14] would be more "QCD-like" in this context and 
it would therefore be interesting to compute the order of the thermal phase transition and 
the Roberge- Weiss endpoint in that model. Recently, it has been demonstrated in QCD 
that the Roberge- Weiss endpoints are triple points for two [42] or three [37] degenerate 
light or heavy flavours, and second order endpoints in the 3d Ising universality class for 
intermediate quark masses. It would be extremely interesting to find a holographic model 
where the triple points turn into second order endpoints, as the quark mass niq is varied. 
This could be realized in models incorporating backreaction of flavours or in defect models 
in general Dp-Dq setups [58]. 
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Appendix A: One-loop free energy at weak coupling 

To calculate the perturbative contributions to the free energy at finite temperature, in the 
presence of an imaginary chemical potential, we need to first understand the perturbative 
effective potential for the Polyakov loop degrees of freedom. The computation is fairly 
standard, so we will only present a brief summary (we refer the reader to [9, 49, 59] for 
details of the theory on fiat space and on S^ x M). At finite temperature, for the SU{N) 
gauge theory, we introduce the Wilson line around the thermal circle as 

U = expi( f Ardr] . (62) 



/o 
Picking a gauge where the zero mode of Ar is diagonal 

A= At dr = Diag(ai, 02, • • • OAf), >^ Oj = mod 27r, (63) 

we want to compute the perturbative contributions to the thermal partition function in 
the presence of these VEVs. Let us suppose that the SU{N) gauge fields are coupled 
to fiA adjoint Weyl fermions and ua real scalars also in the adjoint representation of the 
gauge group. In addition we allow for hf Weyl fermions and nj complex scalars all in the 
fundamental representation of SU{N). The one-loop result for the effective action then has 
the general form 

5i„ioop = ^ logdet (-p2 _ v2(-)) + !^ logdet {-V^ - V^) 



+n/ log det (-Dq - V^ + m^ j - n^ log det {i-f^Vo + i-/'di) 

-hf log det (ij^Do + iYdi - nig) . (64) 



Here, we have allowed for a common mass nig for the fields in the fundamental (/) repre- 
sentation in view of the fact that the theory that we investigate at strong coupling also falls 
in this class. The background field gauge covariant derivatives are defined as 

Vo = dr + iT [A, ] , iDo = dr + iTA, (65) 

when acting on adjoint and fundamental fields respectively. This one-loop effective action 
is then given by 



5i-ioop/Vol (66) 

'■^ dk 



27r2 



k^ 



N 

y (^ (2 + ua) Re In (l - e"'"'^ e''^'') - IfiA Re In f 1 + e"^"'^ e-^^ 



ij=l 



N 

+ Y^ (2nf Re In (l - e"^"' e-'^'^A - 2hf Re In (l + e"^"' e 
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Here 

aij = ai- ttj , Uk = Jk'^ + m^. (67) 

To evaluate the free energy as a function of the imaginary quark number chemical 
potential, we will consider the massless and massive cases separately, as the former can be 
studied exactly at the one-loop order. 

Massless fundamental flavours: With tti = the integrals in Eq. (66) can be done 
to yield the effective potential for the Polyakov loop variables (assuming < «« < 27r and 
^Qj = 0mod27r), 



V^ff 




Oil 
IT 



1 



mod 2 
2" 



(68) 



mod 2 



1 



1 



mod 2 



2n 



/ 



1 



Cti 



vr 



+ 1 



mod 2 



1 



It is easy to see that Qj = Ink/N, k = 0,1,2, ... N — 1 are extrema of the action. The 
contributions from the adjoint degrees of freedom are naturally invariant under Zjv shifts 
Oj — 7- aj + 2iTk/N, while the fundamental flavours break the discrete symmetry. Expanding 
around each of the N extrema, we can show that they are each local minima labelled by 
the integer k, provided 



N{2 + nA + nA) + nf[l 



12^ ) + 2^/ 



1 



6k Qk"^ 



>0. 



(69) 



As long as this condition is satisfied, at the 1-loop level, all N extrema of the effective 
potential will be local minima. While it is interesting to ask what the fate of the effective 
potential is, when the above condition is not satisfied, in the theory that we are interested 
in (namely, AA = 4 SYM coupled to Nf, M = 2 hypermultiplets) - n^ = 6, n^ = 4, 
n^ = rif = 2Nf and N S> Nf ~ the extrema will always be local minima. 

The global minimum of the effective potential is at Ui = 0. 

The effect of the imaginary chemical potential for the fundamental flavours is easily 
captured. Let us, for the sake of simplicity, restrict attention to the case where the complex 
scalars and fermions in the fundamental represention, all have the same (imaginary) chem- 
ical potential. Then, an imaginary chemical potential is introduced by the replacement. 



i?n 



dr + iTA — ifj-i 1 . 



(70) 
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Therefore, the effective potential for /i/ 7^ is the same as (68) with ai — t- ai—l3fij. For small 
non-zero /3/U/, the global minimum continues to be at Oj = 0. However, when 13 m = Inr/N, 
(r G Z) the global minimum shifts to ctj = 2iTr/N. This is simply because the effective 
potential from the fundamental flavours depends on ai — /^/i/, whilst the adjoint sector is 
independent of /U/. One concludes that as /3fij is increased from 0, past n/N, the global 
minimum jumps to Oj = 2it/N following a first order phase transition. The free energy of 
the theory as a function of nj can now be calculated. This is done by first evaluating the 
the effective potential at a given global minimum ai = 0, 

(71) 




NT'^'^ nf 11- { -'-^ -1) I -nf 11 



24 ■' \ V vr mod2 
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This yields the correct free energy, but only in the range — ^ < (3fii < jj- Outside this 
range, for instance when j3^i > ■^, the global minimum moves to Oj = 2tt/N and the free 
energy is given by the function T{[5^i — 2tt/N). The correct free energy as a function of /i/ 
is 

F[/.,] = mm T (/3^, - ^) , (72) 

a function with cusps at /3^/ = (2r — Vjii/N , r G Z. 

When N is large, only the behaviour of F near its local minima is important, i.e., 

T^ / 27rr \ ^ 

F|,4^^^ = (2„, + ,V)iV-,m„f8„-— . (73) 



Massive fundamental flavours: When m ^ 0, the effective potential Eq. (66) cannot 
be expressed in simple analytical form. Nevertheless, it has the same qualitative properties 
as the mq = potential in the presence of non- vanishing fij. Specifically, we can follow the 
reasoning above to deduce the free energy at large A^, 

/bose (y) = Y^ / x"^ dx sinh"^ f -\/x2 + y2 j , (75) 



where 



and 

/fermi iv) = ^ ^^ '^^ COsh"^ f -Vx2+y2j . (Jg) 

As rUq/T — )■ 0, the functions /bose and /fcrmi approach 2 and 1 respectively. 
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